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' Abstract: CMSSM boundary conditions are usually used when calculating cosmological 

dark matter densities. In this paper we calculate the cosmological density of dark matter in 
the MSSM using minimal SOiq soft SUSY breaking boundary conditions. These boundary 
conditions incorporate several attractive features: they are consistent with SOiq Yukawa 
unification, they result in a "natural" inverted scalar mass hierarchy and they reduce 
the dimension 5 operator contribution to the proton decay rate. With regards to dark 
matter, on the other hand, this is to a large extent an unexplored territory with large 
squark and slepton masses mig, large ^0 and small {^,Mx/2}- We find that in most 
regions of parameter space the cosmological density of dark matter is considerably less 
than required by the data. However there is a well-defined, narrow region of parameter 
space which provides the observed relic density of dark matter, as well as a good fit to 
precision electroweak data, including top, bottom and tau masses, and acceptable bounds 
on the branching fraction of Bg fi^ fi^ . We present predictions for Higgs and SUSY 
spectra, the dark matter detection cross section and the branching ratio BR(i35 — > fi^ fi^) 
in this region of parameter space. 
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1. Introduction 

The constrained minimal supersymmetric standard model [CMSSM] |j^ is a well defined 
model for soft SUSY breaking with five independent parameters given by rriQ, Mi Aq, 
tan (3 and sign{^). It has been used extensively for benchmark points for collider searches, 
as well as for astrophysical and dark matter analyses. The economy of parameters in 
this scheme makes it a useful tool for exploring SUSY phenomena. However the CMSSM 
may miss regions of soft SUSY breaking parameter space which give qualitatively different 
predictions. In this paper we consider an alternate scheme, the minimal SO\q supersym- 
metric model [MSOioSM], which is well motivated and opens up a qualitatively new region 
of parameter space. 

In the MSOioSM there are 7 soft SUSY breaking parameters //, Mi/21 ^Oi tan/3, mi% 
(a universal squark and slepton mass), mio (a universal Higgs mass) and Am^ (Higgs 
up/down mass splitting). Moreover the parameters Aq, mio, mi^ must satisfy the con- 
straints [^, m, ^ ~ —2 mi6, mio ^ \/2 mig, m-ig > 1.2 TeV with ji, M^j^ <C m\% and 
tan/3 ~ 50. Note, with these values of the soft SUSY breaking parameters, we can ex- 
plore SUSY phenomena with qualitatively different behavior than in the CMSSM. This is 
mainly due to the Higgs splitting (Am'jj) which, as is well known enables one to obtain 
electroweak symmetry breaking with values mig ^ Also, radiative EWSB with 

tan/3 ~ 50 requires significantly less fine tuning with Higgs mass splitting (see Rattazzi 
and Sarid |^]). Furthermore, with 3 Higgs mass parameters /U, mio, ai^d Am|^ we find 
that the latter two are strongly constrained by EWSB, once we fix the value of fi, which 
we treat as a free parameter. This is unlike the CMSSM where fi is fixed by EWSB. Also 
note that small changes in Am'jj lead to big changes in the CP odd Higgs mass m^ [^. 

It is not at all obvious that the MSOioSM region of soft SUSY breaking parameter 
space is consistent with cosmology Q.^ The dark matter candidate in this model is 

^See also other recent articles discussing Yukawa unification and dark matter M, m 
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the lightest neutrahno. However, since the scalar masses of the first two families are of 
order m-ig > 1.2 TeV, and the third generation sfermions (except for the stops) also tend 
to be heavy, the usually dominant annihilation channels, for the neutralino LSP to light 
fermions via f-channel sfermion exchange, are suppressed. On the other hand, the process 
XX ~^ ff via s-channel A exchange becomes important. This is due to the enhanced CP 
odd Higgs coupling to down-type fermions, which is proportional to tan /?, and because, 
in contrast to heavy scalar exchange, the process is not p-wave suppressed. In an earlier 
analysis, our x^-analysis favored a light CP odd Higgs mass uia ~ 100 GeV 0, although 
heavier A were also allowed. Such light A are however disfavored for two reasons. In order 
to provide efficient annihilation for the LSPs, one would be squeezed into a rather low LSP 
mass region ~ which would require extreme fine-tuning at best. In addition, 

such low rriA are anyway inconsistent with the current limits on BR(i?s ^ I^i this 

analysis, we vary the A mass. 

We study the cosmology of the MSOiqSM in this paper. Obtaining the observed 
relic abundance of cold dark matter, which along with other cosmological parameters has 
recently been determined with an unprecedented accuracy |^] , will provide a new important 
constraint on the model. We also compute the branching ratio for the process Bg ^+ fi^ 
due to A exchange It is absolutely essential to include this latter constraint in our 
analysis. Note, the CDF bound BR{Bs ^ ^+ /U") < 2.6 x 10"^ [|lO|. The cross section for 
the direct detection of dark matter is also computed. In section ^ we define the MSOiqSM, 
describe its virtues and outline the analysis. In section ^ we compute the cosmological 
dark matter density and discuss our results. Then in section Q we discuss our predictions 
for underground dark matter searches and for collider Higgs and SUSY searches. 



2. Minimal SOio SUSY Model - MSOiqSM 
2.1 Framework 

Let us define the minimal SOiq SUSY model. Quarks and leptons of one family reside in 
the 16 dimensional representation, while the two Higgs doublets of the MSSM reside in one 
10 dimensional representation. For the third generation we assume the minimal Yukawa 
coupling term given by 

A 16 10 16. (2.1) 

On the other hand, for the first two generations and for their mixing with the third, we 
assume a hierarchical mass matrix structure due to effective higher dimensional operators. 
Hence the third generation Yukawa couplings satisfy \t = \b = \t = = A. 
Soft SUSY breaking parameters are also consistent with SOiq with 

• a universal gaugino mass M1/2, 

• a universal squark and slepton mass m-ie,^ 

^SOm does not require all sfermions to have the same mass. This however may be enforced by non- 
abelian family symmetries or possibly by the SUSY breaking mechanism. 
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• a universal scalar Higgs mass mio, 

• and a universal A parameter Ag- 
in addition we have the soft SUSY breaking Higgs mass parameters /x and B/j:. Bfi may, as 
in the CMSSM, be exchanged for tan /?. Note, not all of these parameters are independent. 
Indeed, in order to fit the low energy electroweak data, including the third generation 
fermion masses, it has been shown that Aq, ttt-iq, tuiq must satisfy the constraints 

mio « 

mi6 > 1.2 TeV 

/X, Ml/2 < "116 (2.2) 

with 

tan/3ps50. (2.3) 

This result has been confirmed in two recent analyses ||3|, ^.'^ The first property (Eqn. ( |2.2| )) 
is necessary to fit the top, bottom and r masses, in addition to the precision electroweak 
data |2|, 1^, 1^. The second property (Eqn. (|2.3|)) is a consequence of third generation Yukawa 
unification, since mt{mt)/mf,{mt) ~ tan/3. 

One loop threshold corrections at the GUT scale lead to two significant parameters we 
treat as free parameters, although they are calculable in any GUT. The first is a correction 
to gauge coupling unification given by 

eg = [asiMa) - ac] /ac (2.4) 

where the GUT scale Mq is defined as the scale where ai^Mc) = a2{MG) = olq. The 
second is a Higgs splitting mass parameter defined by 

Am^ = (rn\^ - m\^)l2m\^. (2.5) 

In order to fit the low energy data we find eg ~ —4% and Am'jj ~ 13% Q. The largest 
corrections to eg come from the Higgs and SOio breaking sectors, while the correction 
to Amjj is predominantly due to the right-handed r neutrino. For Mp^ ~ 10^^"^"^ GeV 
(appropriate for a light r neutrino mass ~ 0.06 eV) we obtain Amj^ ~ 10 — 7%. 

Finally, as a bonus, these same values of soft SUSY breaking parameters, with mig ^ 
TeV, result in two very interesting consequences. Firstly, it "naturally" produces an in- 
verted scalar mass hierarchy [ISMH] [^]. With an ISMH squarks and sleptons of the first 
two generations obtain mass of order mie at Mz- The stop, sbottom, and stau, on the 
other hand, have mass less than a TeV. An ISMH has two virtues. 

1. It preserves "naturalness" (for values of mig which are not too large), since only the 
third generation squarks and sleptons couple strongly to the Higgs. 

^Note, different regions of parameter space consistent with Yukawa unification have also been discussed 
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2. It ameliorates the SUSY CP and flavor problems, since these constraints on CP 
violating angles or flavor violating squark and slepton masses are strongest for the 
first two generations, yet they are suppressed as 1/mfg. For mig > a few TeV, these 
constraints are weakened IT^ 



Secondly, Super-Kamiokande bounds on r(p — > K^v) > 1.9 x 10 yrs. [14| constrain 



the contribution of dimension 5 baryon and lepton number violating operators. These are 



however minimized with /x, M1/2 *C mie |15 |. 
2.2 Analysis 

We use a top-down approach with a global analysis [16|. The input parameters are 
defined by boundary conditions at the GUT scale. The 11 input parameters at Mq are 
given by — three gauge parameters Mq, aciMc), £3; the Yukawa coupling A, and 7 
soft SUSY breaking parameters fi, M1/2, Aq, tan/3, mfg, mlQ, Arnjj. These are fit 
in a global analysis defined in terms of physical low energy observables. Note we 
keep three parameters {rriiQ, /j,, fixed; while minimizing with the remaining 8 

parameters. Below we will plot contours as a function of fi, M1/2 for different values 
of miQ. We use two (one) loop renormalization group [RG] running for dimensionless 
(dimensionful) parameters from Mq to Mz-'^ We require electroweak symmetry breaking 
using an improved Higgs potential, including and corrections in an effective 2 Higgs 
doublet model below MsusY = \ \{i^'f + "t-? ) 113 • 

The x^ function includes 9 observables; 6 precision electroweak data uem^ Gfj,, ctsi^^z), 
Mz, Mw, Pnew and the 3 fermion masses Mtop-, mb{mi,), M-r- In our analysis we fit the 



central values |T|]: Mz = 91.188 GeV, Mw = 80.419 GeV, G„ x 10^ = 1.1664 GeV 
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'^EAi ~ 137.04, Mj- = 1.7770 GeV with 0.1% numerical uncertainties; and the following 
with the experimental uncertainty in parentheses: as{Mz) = 0.1180 (0.0020), pnew x 10^ = 



-0.200 (1.1) m, Mt = 174.3 (5.1) GeV, mb{mb) = 4.20 (0.20) GeV.^ We include the 
complete one loop threshold corrections at Mz to all observables. In addition we use one 
loop QED and three loop QCD RG running below Mz- 

The output of this analysis is a set of weak scale squark, slepton, gaugino and Higgs 
masses. With regards to the calculated Higgs and sparticle masses, the neutral Higgs 
masses h, H, A are pole masses calculated with the leading top, bottom, stop, sbottom 
loop contributions; while all other sparticle masses are running masses. This output is then 
used to compute the cosmological dark matter density of the lightest neutralino which is 
the LSP. The dark matter analysis is discussed in more detail in section ^. 

Using x^ penalties^ we apply two additional constraints: 



■^Note, we have checked that switching to 2 loop RGEs for dimensionful parameters can be compensated 
for by small changes in the GUT scale parameters, without significant changes in the low energy results. 

^Note we take a conservative error for mtimt) in view of recent claims to much smaller error bars [ p^ . 

^In order to constrain the values of some physical observables in our analysis, such as m^^ or tua, we 
add a significant contribution to the function for values of these observables outside the desired range. 
We refer to this additional contribution as a penalty. Minimization of x^ with Minuit, then pushes the 
fits to the desired range. Of course the x^ penalties then vanish. 
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• > 300 GeV 

• rriA fixed. 

The first is chosen to be consistent with BR(i? ^sj) 0]- Note, although we do 
calculate BR(S Xsj), we do not use it as a constraint in the analysis. This is for two 
reasons — 1) this decay mode depends on 3-2 generation mixing which is model dependent 
and 2) it is not difficult to fit BR(i? — > Xg'y) for values of m^^ > 300 GeV. Hence, in order 
to be generally consistent with the measured value of BR(i? ^sj), we impose m^^ > 300 
GeV. With regards to the second constraint, since and BR(i?5 ^'^ /i^) are both 

sensitive to the value of m^, we fix it's value and present our results for different values of 

Finally, we apply the experimental limits: 

• lower bound on the lightest chargino mass w-^^h- > 104 GeV, 

• lower bound on the light Higgs mass nih > 111 GeV. 

Note, because of the theoretical uncertainty in the calculation of rrih (~ 3 GeV), we conser- 
vatively impose rrih > 111 GeV, instead of the LEP bound for SM Higgs nih > 114.4 GeV. 



3. Cosmological Dark Matter Density 

We compute the relic abundance f^^/i^ of the lightest neutralino using exact expressions 
for neutralino pair annihilation into all allowed final-state channels, which are valid both 
near and further away from resonances and thresholds We further treat the neu- 



tralino coannihilation with the lightest chargino and next-to-lightest neutralino |22] and 



with the lighter stau |£3[ with similar precision. We only neglect the neutralino coanni- 
hilation with the stop which would only affect ^^h? in the regions of parameter space 
which are uninteresting for other reasons, as we comment below. We solve the Boltzmann 

,2 



equation numerically as in [24| and compute with an error of a few per cent, which 



is comparable with today's accuracy on the observational side. The latest determinations 
of cosmological parameters give VlMh? = 0.135j]QQQg for the total matter content and 
Ofe/i^ = 0.0224 lb 0.0009 for the baryonic component. The difference, attributed to cold 
dark matter (CDM), is then 

J7c'£,^^/i2 = 0.113 ±0.009, (3.1) 

which is significantly narrower than previous ranges. We then apply two constraints on 
the dark matter abundance: 

• the upper bound VL^.h'^ < 0.13, 

• 2(7 preferred range 0.095 < Q.yh? < 0.13. 
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mj6=3TeV, m^=500 GeV 




-J r^"^ I I 1 I 1 i I I I I I I I 1 iiiiiiifiifiJ'i 'F™!*^^ 

100 200 300 400 500 600 700 800 



(GeV) 

Figure 1: Contours of constant for ''^w = 3 TeV and tjia — 500 GeV. The red regions are 
excluded by to-^+ < 104 GeV (below and to the left of a black solid curve), mh < 111 GeV (on the 
right) and by fl^h^ > 0.13. To the right of the black broken line one has m/i < 114.4 GeV. The 
green band corresponds to the preferred 2 a range 0.095 < il.^h^ < 0.13, while the white regions 
below it correspond to fl^h^ < 0.095. 



In Figs. 1^-^ we present our results for different values of m-ie and niA in the //, M1/2 
plane. In particular in Fig. |l| we present, for mig = 3 TeV and niA = 500 GeV, the 
(magenta) lines of constant with the cosmologically preferred dark matter region (shaded 
green) satisfying 0.095 < f^^^^ ^ 0.13. We find significant regions of parameter space 
which gives ^ 2, fi^/i^ as above, and satisfies all other phenomenological constraints. 
In addition we have shaded (light red) the regions excluded by collider limits and by 
n^h"^ > 0.13.^ 

In Fig. ^ we present a more detailed analysis of the same mig = 3 TeV, niA = 500 

'^The calculation of BR(_B Xsj) and BR(i3a /x+ fi~) requires a model for fermion mass matrices. 
In the absence of such a model we use the observed CKM matrix elements to calculate these flavor violating 
branching ratios. 

*A similar analysis was performed in the recent paper W. However, they were not able to find acceptable 
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Figure 2: Same as Fig. 1 with contours of constant BK{Bs ^ /i+ /i ) (upper left), nih (upper 
right), n^h^ (lower left) and a^^ (lower right) for mig = 3 TeV and tua — 500 GeV. 



GeV case given in Fig. ||. We now include lines of constant BR(i?s — > /x+ /i^) (upper left), 
ruh (upper right), il^h"^ (lower left), and (lower right). 0"^^ is the spin independent 
neutralino dark matter cross-section relevant for direct dark matter searches. We now 
consider each one of these features further. 

Recall, the branching ratio BR(i?s — > is sensitive to the value of the CP odd 

Higgs mass Q. For rriA = 500 GeV, the branching ratio satisfies 2 x 10"'' < BR(i?s — > 

cosmological solutions. It appears that they did not find any acceptable solutions because they find Yukawa 
unification only for Afi/2 ~ 100 GeV and large ^ ~ 300 GeV. In this region we probably would not find 
acceptable solutions for Q^h^ , no matter what value we take for tua- 
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m,g=2.5TeV, m,=500 GeV m,5=5 TeV, mj=500 GeV 




100 20O 300 40O 600 600 700 100 200 30O 400 500 600 700 800 900 1000 1100 1200 1300 



M,^, (GeV) M,/2 (GeV) 

m,5=2.5TeV, m,= 300 GeV m,g=3 TeV, m^=300 GeV 




M,^, (GeV) M,^, (GeV) 



Figure 3: The same as in Fig. || but for niiQ — 2.5 TeV and niA = 500 GeV (upper left), mig — 5 
TeV and itla = 500 GeV (upper riglit), toiq — 2.5 TeV and niA = 300 GeV (lower left) and mig = 3 
TeV and niA = 300 GeV (lower right). Also marked are contours of constant BR(i?s /i+ /i~). 
The blue regions in the lower two panels are excluded by BR(i?s — > /x+ /i^) > 2.6 x 10^^. Note the 
different mass ranges for M1/2 and fi in the different panels. 



/x+ H^) < 8 X 10^^ for acceptable values of fl^h^ and < 2. In a recent analysis it has 
been shown that, with an integrated luminosity of 15 fb~^, CDF can discover this process 
if BR(i?s jj,^) > 1.2 X 10^^ Hence most acceptable regions of parameter space 

lead to observable rates for BR(i?s fi'^ ^ )• 

In Fig. |2| (upper right) we see that the light Higgs mass increases as M1/2 decreases. 
In the acceptable regions of parameter space we find 116 < nih < 121 GeV. The value of 
the light Higgs mass is however fairly insensitive to mig- 
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The cosmological relic abundance of the neutrahno ^x^"^ (Fig- § (lower left)) is pri- 
marily determined by the direct s-channel pair-annihilation into SM fermion pairs through 
the CP odd Higgs. Since all the sfermions are very heavy, their contribution to reducing 
the neutralino number density is strongly suppressed. In contrast, because of the cou- 
pling Abb (X tan/3 (and similarly for the r's), the ^-resonance is effective and broad. 
Near my, « ruA/^ it reduces il^/i^ down to allowed but uninterestingly small values 
<C 0.1. As one moves away from the resonance, fi^/i^ grows, reaches the preferred range 
0.095 < Q^/i^ < 0.13, before becoming too large f^^^^ ^ 0.13.^ (A similar, but much 
more narrow resonance due to is also present at M1/2 ^ 150 GeV and small fi.) When 
niy. '>^mt {Ml 12 ^420 GeV) and the stops are not too heavy, the LSP pairs annihilate to 
tt-pairs. In the region of large Mi 12, often where is already too low, two additional 
channels become effective. First, in this region the neutralino becomes almost mass degen- 
erate with the lighter stau which leads to reducing Jl^/i^ through coannihilation. Second, if 
ruA is not too large, neutralino pair-annihilation into Higgs boson pairs AA and HH opens 
up. Finally, at /i ^ -^1/2 ) the relic abundance is strongly reduced due to the increasing 
higgsino component of the LSP. 

Finally, the spin independent neutralino cross-section cip^ in the lower right window 
of Fig. m is predominantly determined by the contribution of the heavy CP even scalar 
t-channel exchange to both tree-level and one-loop diagrams. Note that in the preferred 
region of < 2 and 0.095 < n^h^ < 0.13 we find lO'^ pb < a^^ < 10"^ pb. We will 
comment further on our predictions for 0"^^ below. 

In Fig. ^ we display the dependence of our constraints, fi^/i^ and BR(i?s fJ-^ ^ ) 
on rriiQ and tua- By comparing the upper two windows with Fig. |l|we can see that, as mig 
increases, the region with < 2 rapidly grows. Note, the dominant pull in is due to the 
bottom quark mass. In order to fit the data, the total SUSY corrections to mb{mb) must be 
of order —(2 — 4)% In addition there are three dominant contributions to these SUSY 
corrections, a gluino loop contribution oc 03 /i Mg tan/3/m? , a chargino loop contribution 
oc At tan/3/m?^, and a term oc log M|^^y. When miQ increases (with M1/2 fixed) 

the parameter At becomes more negative, since Aq ^ —2 mig and At ~ —3 M1/2 + e ^0 
where e ^ 1. Also, larger values of mig permit a larger range for the ratio m^_^/mi^. Thus 
larger values of mig allows more freedom in parameter space for fitting the data at both 
smaller or larger values of /i, Mx/2- 

In the lower two windows in Fig. |3| we consider two regions with niA = 300 GeV with 
mi6 = 2.5 TeV (lower left) and mie = 3 TeV (lower right). The blue regions are excluded 
by the CDF bound mi{Bs fi+ fi') < 2.6 x 10"^ Note for mie = 2.5 TeV, the region 
with < 2 does not overlap the region with acceptable dark matter abundance (green 
shaded). However for mig = 3 TeV, itla = 300 GeV, Fig. ^ (lower right), we find a small 
region with acceptable x'^ < 1 smd Moreover the branching ratio BR(i?s fi'^ fi^) 

is now close to the CDF bound. On the other hand for mig = 5 TeV, uia = 500 GeV, 
in Fig. ^ (upper right), a new region of parameter space consistent with all the data now 

^Note that at one loop we have Mi(Mz) = M1/2 * ai{Mz)/aG so Mi(Mz) ^ 0.4Mi/2. For bino-Uke 
neutralino (which is true for larger jj,), we thus have ~ 0.4A/i/2- Hence for s-channel annihilation we 
have niA ~ 2m^ « 0.8 M1/2 or M1/2 ~ (5/4) niA for the position of the "peak suppression." 
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opens up with larger /i, This new region becomes cosmologicahy ahowed due to 

neutrahno annihilation into Higgs boson pairs AA and HH and due to coannihilation with 
the lighter stau. 

Hence we see that increasing niA has two effects. It suppresses the branching fraction 
BR(B<j — > /i^). At the same time it moves the s-channel neutralino annihilation 
channel to larger values of M1/2', hence providing larger regions with 0.095 < n^h^ < 0.13 
(compare Figs. || and |3| (lower right) or Figs. ^ (upper and lower left)). In fact, the two 
(green) branches of the preferred range 0.095 < il^/i^ < 0.13 correspond to the two sides 
(except for the upper left window of Fig. |^ where just one side is evident) of the wide A 
resonance in the neutralino pair-annihilation. On the other hand, increasing niA above 
1 TeV or so would move the regions of preferred too far to the right, in potential 

conflict with a lower bound on nih- 

Finally, we comment on BR(i? ^s^)- The current experimental range [26, 27] is 
BR{B Xs7)expt = (3.41 ± 0.36) x 10"^, while the SM prediction, including fuh NLO 
QCD corrections |8|, |9|, is BR{B X^7)sm = (3.70 ± 0.30) x 10"''. In computing the 
SUSY contribution to 6 ^ 57 we further include full LO and dominant NLO-level tan P- 
enhanced contributions Conservatively allowing for the SM+SUSY contribution 

to be in the range (3.41 ± 0.67) x lO"'^ [||, ^ selects a band 300 GeV < < 400 GeV 
which slowly decreases with increasing (This remains approximately true for all the 

cases that we have analyzed except m-ig = 5 TeV and = 500 GeV where one finds a 
narrower range at fi < 200 GeV.) However, BR(i? ^sj) is strongly sensitive to the 2-3 
generation down-type squark mixings [^] which are model dependent and which we do 
not include here. In summary, the process is generally consistent with the most preferred 
regions of M1/2 and n but we do not use it here as a constraint, since it can be easily 
relaxed by employing parameters which are less relevant for our analysis. 

To summarize, we find that regions satisfying all three constraints exist for uiiq > 3 
TeV and > 300 GeV. The acceptable range for fi,Mi/2 grows with increasing rriA 
(for fixed rriiQ = 3 TeV) from approximately 260 GeV < ^ < 295 GeV, 140 GeV < 
Ml/2 < 200 GeV for tua = 300 GeV to 160 GeV < ^ < 520 GeV, 140 GeV < M1/2 < 
530 GeV for uia = 500 GeV. In addition, the allowed regions grow as mie increases 
and for niA = 500 GeV, m-ig = 5 TeV there is also a large fi, M1/2 solution satisfying 
480 GeV < /i < 820 GeV, 640 GeV < M1/2 < 1020 GeV. In the Table we present the input 
and output data from the analysis for three points satisfying all the phenomenological 
constraints. 



4. Predictions and Summary 

In this paper we have analyzed the MSOiqSM and found regions of soft SUSY breaking 
parameter space which fit precision electroweak data, including the top, bottom and tau 
masses and, in addition, fit the cosmological dark matter abundance for the neutralino LSP 
and satisfy BR(i?s — > fj,"^ Z^")- Generically, we find solutions to all the constraints with 
nT'iG > 3 TeV. The squark and slepton masses have an inverted scalar mass hierarchy with 
the first and second generation scalar masses of order mig, while the third generation has 
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mie=3 TeV 



300 400 

(GeV) 



Figure 4: The upper and lower limits on BK{Bs ^ /i+ /i^) as a function of niA in the fi, M1/2 
region of parameter space satisfying all the collider constraints, 0.095 < il^h^ < 0.13 and < 2 
for fixed niie — 3 TeV. 



mass less than 1.3 TeV for mie = 5 TeV. This nice feature of the model suppresses SUSY 
CP and flavor problems. In addition the gaugino masses are typically much lighter, except 
for the large /U, M1/2 region for mig = 5 TeV with a gluino mass of order 1.7 TeV (see 
spectrum in Table for selected acceptable points). 

Note an immediate consequence of such heavy first and second generation sleptons 
is the suppression of the SUSY contribution to the anomalous magnetic moment of the 
muon. We find af^^^ < 2.8 x lO^^'^ (see Table). This is consistent with the most recent 
experimental [^] and theoretical results at la if one uses r-based analysis |34]. However 
it is only consistent with an e^e~-based analysis at 3 a. 

Another interesting result is the enhanced branching ratio for the process Bs — > /^^ 
In Fig. ^ we show the ranges of values of BR(i?s /i^) in the low /x, region of 

parameter space satisfying all the phenomenological constraints with 0.095 < Q.^h? < 0.13 
and < 2 as a function of for fixed mig = 3 TeV. The horizontal red line is the CDF 
bound. Over a significant region of parameter space BR(i?s ^+ //^) > 1 x 10^^ and 
may be observable at the Tevatron (Run II) p5[ |. 

Finally in Fig. ^ we present the cross-section for elastic neutralino-proton scattering 
due to scalar interactions for all regions satisfying the collider constraints, 0.095 < 
< 0.13 and \^ < 2. The green bands are for niiQ = 2.5 TeV, the red for 3 TeV and 
the blue for 5 TeV. The lighter shading is for niA = 300 GeV, the darker for niA = 500 GeV. 
In the last case (mig = 5 TeV and = 500 GeV) there are two branches which correspond 
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Figure 5: Predictions for a^^ vs. for different choices of rriie and niA, subject to the collider 
constraints, 0.095 < fi^/i^ < 0.13 and < 3. The green bands are for mig = 2.5 TeV, the red for 
3 TeV and the blue for 5 TeV. The lighter shading is for tua = 300 GeV, the darker for tjia = 500 
GeV. In the last case (mig — 5 TeV and uia — 500 GeV) there are two branches which correspond 
to the two cosmologically preferred regions in the upper right panel in Fig. |^. 



to the two cosmologically preferred regions in the upper right panel in Fig. |^. Note that 
lower niA generally gives larger a'p^ as expected. For comparison, we also show the bounds 
from the present dark matter searches and the predictions of the general MSSM p5| . 
(Other recent studies of in the case of non-universal Higgs mass in a variant of the 
CMSSM can be found in |^6[.) Over the next two to five years the experimental sensitivity 
is expected to gradually improve by some three orders of magnitude. This will cover large 
parts of the predicted ranges of dp^, especially at lower values of niA where Bs — > /i~ 
will be accessible at the Tevatron (Run II). 
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